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The study of electron transport and scattering processes limiting electron mobility in high-quality 
semiconductor structures is central to solid-state electronics. Here, we uncover an unavoidable 
source of electron scattering which is caused by fluctuations of nuclear spins. We calculate the 
momentum relaxation time of electrons in quantum wells governed by the hyperhne interaction 
between electrons and nuclei and show that this time drastically depends on the spatial correlation 
of nuclear spins. Moreover, the scattering processes accompanied by a spin flip are a source of the 
backscattering of Dirac fermions at conducting surfaces of topological insulators. 
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I. INTRODUCTION 


The invention of modulation-doped semiconductor 
structures [I| and the subsequent progress in semiconduc¬ 
tor technology have led to the fabrication of ultrahigh- 
mobility two-dimensional electron systems ii and dis¬ 
covery of novel exciting quantum phenomena such as the 
fractional quantum Hall effect and microwave-induced 
resistance oscillations Much effort is focused now on 
the search for new technological approaches and the op¬ 
timization of quantum well (QW) design to reduce struc¬ 
ture disorder and increase electron mobility. This raises 
the question of the fundamental limitation of electron 
mobility that could be achieved in defect-free QW struc¬ 
tures with ideal interfaces, see Ref. 0 for a recent dis¬ 
cussion. 

Here, we analyze a source of electron scattering stem¬ 
ming from hyperfine interaction between electron spins 
and spins of nuclei constituting the crystal lattice. This 
scattering mechanism is unavoidable in HI-V compounds 
since all stable and long-lived isotopes of anions (N, P, As, 
Sb) and cations (B, Al, Ga, In) possess non-zero nuclear 
spins. Hyperfine interaction in semiconductors has been 
extensively studied in the context of coupled electron 
and nuclear spin dynamics in bulk semiconductors Q, 
quantum dots (QDs) @,1^ and QWs [lol - [l^ . and also in 
the spin-dependent electron transport along edge chan¬ 
nels of a two-dimensional electron gas through 

QDs (spin-blockade effect) [20l - l^ , but not in bulk charge 
transport measurements. Here, we calculate the electron 
mobility limited by electron-nuclear hyperfine interaction 
in QWs for various spin configurations of the electron and 
nuclear subsystems including the case of unpolarized elec¬ 
trons and nuclei and the dynamic nuclear polarization 
(DNP). It is shown that the effect of electron-nuclear in¬ 
teraction on the electron mobility drastically depends on 
the spatial correlation of nuclear spins. Generally, both 
spin-conserving and spin-flip processes contribute to the 
electron scattering. The quadrupole splitting of the nu¬ 


clear spin levels in strained QWs or the Zeeman split¬ 
ting of the electron and nuclear levels in a magnetic field 
can suppress the spin-flip electron scattering by nuclear 
spin fluctuations at low temperatures. A uniform nuclear 
polarization achieved, e.g., by DNP, suppresses in turn 
the spin-conserving scattering processes. Spin-flip scat¬ 
tering is an unavoidable source of the backscattering of 
two-dimensional Dirac fermions emerging at conducting 
surfaces of topological insulators. 


II. SCATTERING BY UNCORRELATED 
NUCLEAR SPINS 


The effective Hamiltonian of hyperfine interaction be¬ 
tween the conduction-band electrons and nuclei can be 
presented in the form Q 

V = '^AaVo S ■ Ia,nS{r - Ra,n) , (1) 

a.,n 


where a is the index of nucleus species, n enumerates 
nuclei of certain species, Aa are the constants of interac¬ 
tion, vq is the volume of the primitive cell, S and la^n are 
the electron and nucleus spin operators, respectivel y, a nd 
Ra,n = {Pa,n, Za,n) are the positions of the nuclei |24l |. 

We assume first that nuclei are unpolarized on average 
and their spin states are uncorrelated with each other. 
Then, the scattering of an electron by different nuclei oc¬ 
curs independently and the total probability of the scat¬ 
tering from the initial state (fc, s) to the final state (fc', s') 
is given by the golden rule rate 


X ^ {Sk + Es + Eaj — 6k' — 6s' — Eaj') , 


(a,n) |2 


sj I 


Paj 


( 2 ) 


where k and k' are the wave vectors in the QW plane, 
s,s' = ±1/2 are the electron spin projections, j and j' 
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are the initial and final nuclear spin projections, 
is the matrix element of scattering at the potential given 
in Eq. Paj is the nuclear spin distribution function, 
£k = is the electron kinetic energy, m* is 

the effective mass. Eg is the energy related to electron 
spin, e.g., in an external magnetic field, and Saj are the 
nuclear energies. The hyperfine interaction is weak and 
short-range as compared to the de Broglie wavelength of 
electrons. It allows us to study the scattering in the first 
Born approximation and neglect screening. 

The squared modulus of the scattering matrix element 
has the form 

= ^1^0 l(s'/|5' • Ia,n\sj)\^ , (3) 

where tpi^) is the function of electron size quantization 
in the QW. The operator S ■ I can be rewritten in the 
form 

S-I = S,h + iS+I-+S-I+)/2, (4) 

where S± = Sx ^ iSy, I± = /a, ± i/y, and S 13 and Ijs 
(/3 = X, y, z) are the Cartesian components, which yields 

{s'f\S-I\sj) =sj5s's5j'j 
+ \\/ {I ~ j)il + J + 1) Ss' ,s-lSj' J+l 
+ - j + 1)(.^ + 3)^s',s+i5j',j-i ■ (5) 

The electron relaxation time Tp determining the mobil¬ 
ity is expressed via the probability of scattering ([2). To 
derive the expression for Tp we follow Boltzmann’s ap¬ 
proach and write down the collision integral, which plays 
the role of the “friction force” in Boltzmann’s equation, 

St/fcs = 'y^\fk's’{^ — fks)Wks,k's' — fks{^ — fk's')Wk’s',ks\^ 

k' s' 

(6) 

where fks is the electron distribution function. In the 
presence of a weak driving electric field, the distribution 
function has the form fks = fs^\£k) + Sfks, where /i°^ is 
the equilibrium function and Sfks is a small anisotropic 
correction. To first order in Sfks and for Wks,k's' inde¬ 
pendent of the directions of the wave vectors k and fc', 
which is valid for short-range scattering by nuclei, Eq. (|6]) 
yields 

Stfks = -^, ( 7 ) 

'^P 

where 

=Y.^fi'\sk')Wks,k's' + [l-f^?\£k')]Wk's',ks}- ( 8 ) 

k' s' 

Finally, for the scattering probability given by Eq. m, 
homogeneously distributed nuclei and degenerate elec¬ 
tron gas, we obtain 

\{s'j'\S ■ Ia\sj)\^ 

a jj's' 

^ [Paj {Paj' Poi-j^fs' i.^F T £otj £s' £aj'')] 

X Oi^Ep £aj £s' £aj'^ ; ( 9 ) 


where IVa are the densities of nuclei of certain species, 
Ep is the Fermi energy, and 9(e) is the Heaviside step 
function. We note that, for rectangular QWs with in¬ 
finitely high barriers, J ip'^{z)dz = 3/{2d) with d being 
the QW width. Also note that, in general, Tp can depend 
on the electron spin s unless the nuclear-spin distribution 
is symmetric, pj = P-j. The time reversal symmetry 
imposing the condition pj = p-j can be broken by an 
external magnetic field or DNP (studied below), spon¬ 
taneously (see discussions on the possibility of nuclear 
self-polarization in Refs. Bill), or, in principle, by driv¬ 
ing an electric current through the quantum well. In the 
latter case, the electric current induces a spin polariza¬ 
tion of electrons due to spin-orbit coupling which results, 
in turn, in a build up of DNP [ 2 ^. The emerging elec¬ 
tron and nuclear spin polarizations are both proportional 
to the electric current and, therefore, do not affect the 
linear electron transport. Non-linear effects are beyond 
the scope of this paper. Below, we analyze Eq. (I9|) for 
some particular cases which can be of interest. 

(i) All spin states of nuclei are degenerate and equally 
populated, £aj = 0, Paj = ll( 2 Ia + 1), where are 
the nuclear spins; electrons are unpolarized, £s = 0. In 
this case, both spin-conserving and spin-flip scattering 
processes are allowed. Taking into account that 

^Ksj|5-J|sj)p=J(/+l)(2/+l)/12, 

j 

^ Ks ± 1 , J T 1 |^ • I\s,j)\^ = /(/ + 1 )( 2 / + l)/6 , 

j 

we obtain 

a. 

We note that spin-flip scattering processes in the context 
of electron and nuclear spin relaxation were theoretically 
considered in Ref. d 

The degeneracy of nuclear spin levels may be lifted due 
to the quadrupole interaction of the nuclear spins with 
the strain-induced gradient of the crystal field in lattice- 
mismatched structures Q. In a simple axial model rele¬ 
vant for (OOl)-oriented QWs, the quadrupole interaction 
determining the splitting and order of the spin levels is 
proportional to i.e., levels with different p have dif¬ 
ferent energies. The sign of the quadrupole splitting is 
opposite for axial tension and compression. Typical val¬ 
ues of the quadrupole splitting in III-V heterostructures 
are of the order of 1 - 10 neV which corresponds to 10“® - 
10“^ K at the temperature scale [1, [l^. Therefore, non¬ 
equal thermal population of the nuclear spin levels at 
sub-mK temperatures can occur and affect the electron- 
nuclear interaction. 

(ii) Nuclear spin levels are split by strain in such a 
way that the ground levels are characterized by the high¬ 
est spin projections j = ±/a, la > 1/2; electrons are 
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unpolarized, = 0. Temperature is lower than the nu¬ 
clear quadrupole splitting and, therefore, only ground 
levels contribute to scattering, paj = 1/2 if j = ±Jq, 
and Paj = 0 otherwise. Only spin-conserving scatter¬ 
ing processes can occur and the corresponding electron 
relaxation time is given by 

r-^ = S / Nall. (11) 

(iii) Nuclear spin levels are split by strain in such a 
way that the ground levels are characterized by the low¬ 
est spin projections j = ±1/2, la is half-integer; electrons 
are unpolarized. Eg = 0. Temperature is lower than the 
quadrupole splitting energy and only ground levels con¬ 
tribute to scattering, Pa,±ij 2 = 1/2- In this case, both 
spin-flip and spin-conserving processes contribute to scat¬ 
tering and the electron relaxation time has the form 

T-p^ = ^ /Alvl Na[l/2 + {la + 1/2)^]. 

( 12 ) 

(iv) All spin states of nuclei are degenerate and equally 
populated, Eaj = 0, Paj = 1/(2/q ± 1); electrons are 
completely spin polarized by an external magnetic field 
(thermal nuclear polarization is small and neglected). In 
this particular case, only spin-conserving processes can 
occur and Eq. (|9|) yields 

/ ^*iz)dz'^AlvlNaIa{Ia + l) ■ (13) 


In 111-V semiconductor structures, nuclear spins can be 
efficiently polarized by DNP [ 3,0 ■ In the case of uniform 
nuclear polarization, the average nuclear field should be 
excluded from the Hamiltonian of electron-nuclear in¬ 
teraction causing the electron scattering because it does 
not introduce any disorder which breaks the translational 
symmetry of the crystal lattice. Accordingly, the opera¬ 
tors la in Eqs. m and (|ni) should be replaced by la — la, 
where the vectors la are given by 


f _ Y.I3^0NpIp 

““ A.Y.pNp 


(14) 


the index /3 runs over the nuclei of anions or cations if a 
stands for a nucleus of an anion or a cation, respectively, 
and la is the average nuclear spin of a certain isotope. 
In the simple case of an isotopically pure crystal, where 
all anions and cations are of certain isotopes, la = la- 
We assume that all the vectors la point along the same 
axis 2 and choose z as the spin quantization axis. Then, 
one obtains la = ja = J2j jPaj and the scattering marix 
elements 


{s'j'\S ■ {la - ia)\sj) =s{j- ia) Sg'sSfj 

+ 5'\/ {la — j){Ia ± J ± 1) 5s' ,s-ldj' J+1 

+ ^\/{Ia — j + ^){Ia + j) 5s',s+l5j'J-l . (15) 


Now we discuss the electron relaxation time in the pres¬ 
ence of DNP. 

(v) The nuclei are spin polarized; electrons of both 
spin states are present at the Eermi level. The splitting 
of nuclear levels is lower than temperature so that both 
spin-conserving and spin-flip processes can contribute to 
scattering. In such conditions, the momentum relaxation 
time becomes spin dependent and is given by 

rp,l = ^J'll^^{z)dz^ AlvlNa[ll+Ia+Ia-‘ 2 ja{L + s)]- 

Ot 

( 16 ) 

(vi) The nuclei are spin polarized; electrons are com¬ 
pletely spin polarized by the nuclear (Overhauser) effec¬ 
tive magnetic field or an external magnetic field. Then, 
only spin-conserving scattering occur and we obtain the 
relaxation time 

Tp^ = ^ Alvl Na[jl - 2jaia + /J] , 

_ (17) 

where ja = Jljj^Paj- For isotopically purified crys¬ 
tals, the expression in the square brackets is reduced to 
ja — Ja • This quantity is zero for fully polarized nuclei 
when both and j^ are equal to I^. The absence of 
spin-conserving electron scattering by the fully polarized 
nuclei, when the crystal translational symmetry is re¬ 
stored, is in accordance with the free motion of Bloch 
electrons in a periodic potential. In isotopically mixed 
QW structures, the scattering may occur even in the case 
of fully polarized nuclei due to the disorder caused by a 
difference in the nuclear spins and/or the interaction 
constants Aa of anions (or cations). The corresponding 
relaxation time is given by 

= ^/ jj^{z)dzJ2 Alvl Na{Ia - Lf ■ ( 18 ) 

CK 

Now we estimate the electron relaxation time gov¬ 
erned by the scattering from unpolarized nuclei follow¬ 
ing Eq. (H®. The estimation for a 10-nm-wide QW 
grown from GaAs, where m* ~ O.OGTttiq with mo be¬ 
ing the free electron mass, Ica = 3/2, /as = 3/2, 
and vlJ^aA^Na ~ 0.2 x 10“^'^ meV^ cm^ (Ref. 1^ . 
gives Tp ~ 10“® s. This corresponds to the mobility 
/i ^ 4 X 10^^ cm^/(Vs) that is still a few orders of magni¬ 
tude higher than the mobility achieved to date 0. The 
relaxation time governed by the hyperfine interaction can 
be much shorter in structures made of atoms with large 
nuclear spins (e.g., An = 9/2) or heavy atoms, where the 
interaction constants are larger. 

The probabilities of spin-conserving and spin-flip pro¬ 
cesses are comparable. Therefore, electron scattering by 
unpolarized nuclei makes a contribution to electron spin 
relaxation with the time Tg comparable to the momen¬ 
tum relaxation time Tp calculated above. This spin re¬ 
laxation mechanism can be important if other mecha¬ 
nisms are suppressed, e.g., in (llO)-oriented QWs where 
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the D’yakonov-Perel’ spin relaxation mechanisms is sup¬ 
pressed for the out-of-plane spin component [2§-[3lj and 
the spin lifetime up to 0.5 /iS has been recently deter¬ 
mined [111. 


where (pk is the polar angle of the vector k. Following 
Eqs. HI) and (O one can readily calculate the matrix 
element of scattering. For the scattering from the initial 
state k to the final state k' the squared modulus of the 
scattering matrix element has the form 


III. SCATTERING OF DIRAC FERMIONS 


Spin-flip processes contribute to the backscattering of 
two-dimensional Dirac fermions emerging at surfaces of 
3D topological insulators (TIs). In such materials, the 
strong spin-orbit interaction giving rise to topologically 
protected surface states locks the carrier spin and mo¬ 
mentum [s^. Therefore, elastic scattering between the 
states with the opposite momenta is forbidden in the 
presence of time reversal symmetry. Interaction with nu¬ 
clear spins breaks the time reversal symmetry in the sub¬ 
system of the Dirac fermions and enables the backscat¬ 
tering. Electron-nuclear interaction leading to backscat¬ 
tering between one-dimensional helical edge states of 2D 
TIs was studied in Refs. in and [li Here, we calculate 
the probability of scattering for two-dimensional Dirac 
fermions on surfaces of 3D TIs and provide estimates for 
HgTe and Bi 2 Se 3 . 

HgTe is a gapless semiconductor with the inverted 
band structure which becomes a 3D TI if strained and 
the strain opens a gap in the otherwise four-fold degen¬ 
erate Fg states [^.1^. Within the 6-band k-p theory, 
relevant for narrow-band materials, the topological sur¬ 
face states are described by the wave functions 




m=±l/2 


(19) 


m=±l/2,±3/2 


exp(ife • p ), 


where ipTe,m{z) and 'ipTs,miz) are the envelope functions 
in the direction normal to the surface, iFg, m) and jFg, m) 
are the the basis Bloch amplitudes of the s-type Fg and 
p-type Fg states, respectively, and k is the wave vector 
in the surface plane. For strained HgTe films, the wave 
function m contains considerable contribution (~20 %) 
of the Fg states [S^. Since the hyperfine interaction 
for s-type Bloch amplitudes given by the Fermi contact 
term is much stronger than that for p-type Bloch am¬ 
plitudes [1, , a good estimation is that the electron- 

nuclear interaction for the Dirac fermions at HgTe surface 
is of contact type and determined by the contribution of 
the Fg states. 

The electron spin of the surface states in the axial ap¬ 
proximation lies in the surface plane and points perpen¬ 
dicular to the wave vector. Accordingly, the spinor com¬ 
posed of the functions 4 ’re,±i/ 2 iz) can be presented in 
the form 


i^re, + l/2{z) 

1 

1 


“72 

iexp{ipk) 




( 20 ) 


+ {Ia — j + l)(da + + {la — j){Ia + j + l)<^i'J-l-l], 

where 9 = pk' ~ Vk is the angle of scattering. 

The differential probability of elastic scattering of a 
test particle by the angle 9 in 2D systems is given by 

^ 2ttWv ^ ^ 

where v = {l/h)dek/dk is the velocity which is indepen¬ 
dent of the energy for linearly dispersive Dirac fermions. 

We assume that all spin states of the nuclei are degen¬ 
erate and equally populated, Paj = 1/(2/q, -|- 1). Then, 
summing up over the nuclei, we obtain 

= 24^ + sin2(0/2)] d9 

X [ \i^rA^)\U^Y.^^o.^oNaIa{Ia + l)- (22) 


Equation describes the scattering of two-dimensional 
Dirac fermions by nuclear spin fluctuations. It shows 
that the backscattering {9 = tt) determined by spin-flip 
processes is twice as efficient as the forward scattering 
{9 = 0) determined by spin-conserving processes. 

In natural HgTe, about 30 % of H g nu clei and about 
8 % of Te nuclei possess non-zero spins [l^ • To the best of 
our knowledge, the hyperfine interaction constants have 
not been measured yet. Considering the fraction of nuclei 
with non-zero spins, the typical parameters of surface 
states in HgTe films [s^: k = 2 x 10® cm“^, v = 0.5 x 
10® cm/s, the characteristic length of the surface-state 
localization d = 10 nm, and the Fg band partition 0.2, 
we estimate that the probability of scattering by nuclear 
spin fluctuations is two-three orders of magnitude lower 
than that in GaAs quantum wells. 

Other prominent examples of 3D TIs are binary and 
ternary compounds of Bi with Se and Te [s^. Natural Bi 
consists of the only isotope ^®®Bi with the nuclear spin 
9/2. Recent measurements of nuclear magnetic resonance 
in n-type Bi 2 Se 3 crystals have revealed very strong con¬ 
tact interaction between electrons and ^®®Bi nuclei (in 
spite of the fact that the Bloch amplitude is mostly of 
p-type) [ 33 . The contact hyperfine interaction constant 
is found to be comparable and even exceed those for the 
conduction-band electrons in GaAs. These results sug¬ 
gest that the scattering of Dirac fermions at the Bi 2 Se 3 
surface by nuclear spin fluctuations can be as efficient as 
that for electrons in GaAs-based structures. 
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IV. SCATTERING BY FLUCTUATIONS OF 

MACROSCOPIC NUCLEAR POLARIZATION 

The electron-nuclear interaction is drastically en¬ 
hanced if the nuclear spins are spatially correlated and 
polarized at a macroscopic scale. Such a nuclear spin 
polarization inhomogeneous in the QW plane can be cre¬ 
ated, e.g., via DNP by optical grating technique [3§ - l40l| . 
In the case of macroscopic nuclear polarization, the in¬ 
teraction can be described as the Zeeman term 

V{p) = gpBS- B^ip) (23) 

with the effective nuclear (Overhauser) magnetic field 

Bn{p) = ^{AaVoNa/gPB) f Ia{P, z) {z)dz (24) 

cx. 

which varies in the QW plane at a scale much larger 
than the crystal lattice constant. Here, is the Bohr 
magneton, g is the effective electron g-factor, and the 
overline denotes quantum mechanical averaging over the 
ensemble of nuclear wave functions. 

The effective magnetic field Eq. (1^ produces a spin- 
dependent electron potential which causes the scattering. 
For the effective field Bn{p) oriented along a certain axis, 
e.g., the growth direction, the momentum relaxation time 
assumes the form 

^ ^ y (J5n(p)-Bn(p'))2fef.sine/2(l-COs6»)d6», 

(25) 

where {Bn{p) ■ JB„(p'))q is the Fourier image of the spa¬ 
tial correlation function (JB„(p)-JB„(p')), kp is the Fermi 
wave vector, and 9 is the angle of scattering. 


To estimate the momentum relaxation time we assume 
that the nuclear polarization is randomly distributed in 
the QW plane with the characteristic correlation length I 
and zero mean value. Then, for GaAs-based QWs with 
the nuclear polarization la ~ 1, vq ^aNa ~ 0.1 meV 
(Ref. [ 2 ^, the Fermi wave vector kp = 10® cm“^, and the 
correlation length I ^ 1/kp, one obtains Tp ^ 10“® s. 
Such Tp is comparable to the momentum relaxation time 
in high-mobility structures. It indicates that strong spa¬ 
tially inhomogeneous spin polarization of nuclei can con¬ 
siderably affect the electron transport in quantum wells 
and can be probed by electrical measurements. Spatially 
oscillating nuclear polarization created by optical grating 
technique may cause the Bragg diffraction of electrons 
and has even stronger impact on the electron transport. 

To summarize, we have calculated the limitation of 
electron mobility in III-V quantum wells from the un¬ 
avoidable source of disorder originating from the fluc¬ 
tuations of nuclear spins. We have analyzed various 
spin configurations of the electron and nuclear subsys¬ 
tems and shown that the electron mobility determined 
by the electron-nuclear hyperfine interaction drastically 
depends on the spatial correlation of nuclear spins. While 
the electron mobility limited by the hyperfine interaction 
with uncorrelated nuclear spins is still few orders of mag¬ 
nitude higher than that achieved in high-mobility quan¬ 
tum wells, nuclear spins that are spatially correlated and 
polarized at a macroscopic scale can considerably affect 
the electron transport in modern high-mobility quantum 
wells. 
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